Influence of electron-acoustic phonon scattering on intensity power broadening in a 

coherently driven quantum-dot cavity system 
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We present a quantum optics formalism to study intensity power broadening of a semiconductor 
quantum dot interacting with an acoustic phonon bath and a high Q microcavity. Power broadening 
is investigated using a time-convolutionless master equation in the polaron frame which allows for 
a nonperturbative treatment of the interaction of the quantum dot with the phonon reservoir. We 
calculate the full non-Lorentzian photoluminescence (PL) lineshapes and numerically extract the 
intensity linewidths of the quantum dot exciton and the cavity mode as a function of pump rate and 
temperature. For increasing field strengths, multiphonon and multiphoton effects are found to be 
important, even for phonon bath temperatures as low as 4 K. We show that the interaction of the 
quantum dot with the phonon reservoir introduces pronounced features in the power broadened PL 
lineshape, enabling one to observe clear signatures of electron-phonon scattering. The PL lineshapes 
from cavity pumping and exciton pumping are found to be distinctly different, primarily since the 
latter is excited through the exciton-phonon reservoir. To help explain the underlying physics of 
phonon scattering on the power broadened lineshape, an effective phonon Lindblad master equation 
derived from the full time-convolutionless master equation is introduced; we identify and calculate 
distinct Lindblad scattering contributions from electron-phonon interactions, including effects such 
as excitation-induced dephasing, incoherent exciton excitation and exciton- cavity feeding. Our 
effective phonon master equation is shown to reproduce the full intensity PL and the phonon- 
coupling effects very well, suggesting that its general Lindblad form may find widespread use in 
semiconductor cavity-QED. 

PACS numbers: 42.50.Ct, 78.67.Hc, 78.55.-m 



I. INTRODUCTION 

Semiconductor quantum dots (QDs) embedded in 
microcavities have established themselves as a new 
paradigm in cavity quantum electrodynamics (cavity- 
QED). Technological progress in the design and fabrica- 
tion of semiconductor cavity-QED systems has enabled 
them to be used as components in quantum information 
processing^^^ and for the generation of indistinguishable 
photonS^^. These quantum applications require robust 
cavity-QED based QD devices that rest on the ability 
to manipulate and control the underlying quantum pro- 
cesses. Such quantum control is usually obtained when 
the cavity and QD are in the intermediate to strong cou- 
pling regime^— . Recent experimental studies have fo- 
cussed on resonance fluorescence of a QD coupled to a 
cavity modeiSrJ^, and significant progress has been made 
in the study of an off-resonant QD cavity system that is 
used to observe resonance fluorescence of a single pho- 
ton emitter— ii^. Semiconductor micropillar systems are 
particularly attractive since a geometrical separation be- 
tween the pump field and emitted fluorescence signal can 
be made, facilitating nonlinear quantum optical studies 
such as intensity power broadeningi^. 

For semiconductor cavity-QED systems, signatures of 
acoustic phonon processes have been noted with inco- 
herent excitation, resulting in off-resonant "cavity feed- 
•^g57i5 22 asymmetric {on-resonance) vacuum 

Rabi double t^^i^^i^^ . Various phonon-coupling mod- 
els have been developed to try and explain these fea- 



ture a^^i^^ — ; for example, data obtained for the linear 
spectrum of single site-selected dots in cavities show good 
agreement with photon Green function theories — where 
the phonon coupling is included as a self-energy cor- 
rection to the spectrun>22ri24. Recently, several works 
have also experimentally investigated coherent power (in- 
tensity) broadening in semiconductor cavity-QED sys- 
tems. For example, Majumdar et. al^^ studied the role 
of phonon-mediated dot-cavity coupling on the power 
broadened PL intensity for a planar photonic crystal sys- 
tem; with experiments performed at temperatures of 30- 
55 K on self-assembled InAs QDs, the cavity-emitted in- 
tensity PL were found to have extraneously broadened 
linewidths relative to a bare QD (when compared with 
calculations from a simple atomic ME). While additional 
coupling may occur, e.g., from the QD to the contin- 
uum states due to the presence of the nearby wetting 
layer (if it exists), or due to Auger scattering^^, these 
processes are usually more important for incoherent ex- 
citation. For near-resonant coherent excitation, Ulhaq 
et. al^ have demonstrated that dephasing and coupling 
due to acoustic phonons is likely the primary (and in- 
trinsic) mechanism that couples the QD and the spec- 
trally detuned cavity mode; their experiments were per- 
formed using self- assembled InGaAs/GaAs QDs embed- 
ded in a single cavity layer of a semiconductor micropil- 
lar cavity. While the important role of electron-phonon 
scattering on the linear absorption and emission spec- 
tra of self-assembled QDs is now becoming better estab- 
lished^—, there appears to be little theoretical work 
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describing phonon effects on PL power broadening in a 
semiconductor cavity-QED system. 

Nonlinear resonance fluorescence of an InGaAs QD 
embedded in a high-quality micropillar cavity was re- 
cently investigated by Ulrich et al.^^ ^ where, in contrast 
to atomic cavity-QED, a clear indication of excitation- 
induced dephasing (EID) was found to manifest in a Mol- 
low triplet spectra with pump-induced spectral sideband 
broadening. Without cavity interactions, it is well known 
that the interaction of the driven QD with the underlying 
phonon reservoir can introduce additional dephasing pro- 
cesses and acoustic phonon sideband^—. In a cavity 
system, these phonon processes can also result in signif- 
icant coupling between a non-resonant cavity and a QD 
exciton. Very recently, a polaron master equation (ME) 
description of phonon-induced EID in QDs and cavity- 
QED was described by Roy and Hughes^; McCutcheon 
and Nazir have also adopted a polaron ME approach 
to describe pule-excited excitons (without cavity inter- 
actions)^. In light of these phonon scattering studies 
and the emerging class of semiconductor cavity-QED ex- 
periments, the inclusion of phonon-scattering in the the- 
oretical description of PL power broadening in a cavity- 
QD system is highly desired. More generally, one desires 
accurate quantum optical descriptions of the semicon- 
ductor cavity-QED system, where important electron- 
phonon interactions are accounted for. 

In this paper, we present a quantum ME formalism 
to study the intensity power broadening of a semicon- 
ductor cavity-QED system and identify the qualitative 
features of power broadening in the intensity PL intro- 
duced due to electron-phonon interactions. We exploit a 
time-convolutionless ME (i.e., local in time) for the re- 
duced density matrix of the dot-cavity subsystem, where 
the system-bath incoherent interaction is treated to sec- 
ond order— The perturbative treatment is performed 
in the polaron frame which allows us to study the ef- 
fects of phonon dephasing on the coherent part of the 
Hamiltonian exactly. Importantly, the cw laser driving 
the QD introduces additional EID effects in addition to 
pure dephasing due to the phonon reservoi r "^^i^^ . In the 
appropriate limits, the model fully recovers the indepen- 
dent boson model (IBM)^ — and the Jaynes-Cummings 
model. The polaron transform is particularly convenient 
for studying QD cavity-QED systems as it eliminates the 
exciton-phonon coupling and introduces a modified dot- 
cavity coupling and a modified radiative decay rate^J^^; 
in addition, there is a phonon-induced renormalization of 
the QD resonance frequency through the polaron shift. In 
the case of an exciton driven system, the Rabi frequency 
of the cw laser is also renormalized by a temperature- 
dependent factor which essentially accounts for the de- 
phasing of the cw drive due to phonon coupling. A similar 
polaron ME approach was previously derived by Wilson- 
Rae and Imamoglii^^, who studied the linear absorption 
spectrum of a cavity-QED system; however, their ME 
form^ was non-local in time and is substantially more 
difficult to solve than the time-convolutionless forn>^. 



For our coherently-pumped cavity-QED investigations, 
we consider two distinctly different pumping scenarios: 
(i) the QD is driven by a coherent continuous wave (cw) 
laser field, and (n) the cavity mode is driven by a coher- 
ent cw laser field. We describe the generic features arising 
due to the relative interplay between phonon-induced dot 
cavity coupling and EID in the case of a QD driven sys- 
tem, and compare and contrast with power broadening 
for a cavity driven system; we also discuss the differences 
in the integrated PL (IPL) for a dot-driven and cavity- 
driven system. For strong coherent drives (fields), the 
intensity PL contain significant phonon bath signatures 
over a wide range of frequencies. To help explain the ef- 
fects of phonon scattering in these systems, we also derive 
an effective phonon ME, of the Lindblad form which is 
shown, in certain regimes, to yield very good agreement 
with the full polaron ME. 

Our paper is organized as follows. In SecHIlwe present 
the model Hamiltonian and derive a time-convolutionless 
polaron ME where electron-phonon interactions are in- 
cluded to all orders. In Sec. IIIII we introduce an effec- 
tive phonon- modified Lindblad ME and compare it to 
the full time-convolutionless solution; the effective Lind- 
blad ME is shown to yield good agreement with the time- 
convolutionless ME, and we use it to describe the various 
phonon scattering processes. In Sec. llVf A-E) we present 
and discuss our numerical results of the power broadening 
lineshape for both QD-driven and cavity-driven systems. 
In Sec. Owe present our conclusions. Appendix [X] pro- 
vides some technical details about the derivation of our 
effective phonon scattering rates and Lindblad ME. 



II. GENERAL THEORY AND POLARON 
MASTER EQUATION MODEL 

The dynamics of a strongly confined QD can be mod- 
eled by considering a quantized electron-hole excitation, 
where the electron occupies a conduction band state and 
the hole occupies a valence band state. Neglecting quan- 
tum spin, the dominant features of a strongly confined 
QD can be described by the two lowest energy bound 
states. This two- level model is then conditioned by the 
interaction of the electrons with the lattice modes of 
vibration, i.e., the acoustic phonons. When the effec- 
tive two-level system is driven by a cw laser field, power 
broadening may be substantially modified by the cou- 
pling of the QD to the phonon modes^i^. Figure [1] 
shows a schematic of a semiconductor cavity-QED sys- 
tem [Fig. HJa)], and an energy- level diagram associated 
with cavity-pumping [Fig. mb)] and exciton pumping 
[Fig. IHc)]; the various parameters in the figures will be 
introduced below. The semiconductor cavity system of 
interest could be a micropillar cavity system [cf. Fig.[Ha)] 
which allows one to excite and measure through differ- 
ent photon reservoirs^S (e.g., cavity pumping and exciton 
emission). 

Working in a frame rotating with respect to the laser 
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FIG. 1: (a) (Color online) Schematic of an example semi- 
conductor cavity used in cavity-QED (micropillar system), 
containing a coupled QD and driven by a cw laser {rfx^ from 
the side, or ?7c, from the top). The micropillar system is ad- 
vantageous to study PL power broadening since, e.g., exciton 
broadening can be directly studied through cavity emission — 
since these input /output channels are geometrically decou- 
pled, (b) Schematic of a cavity driven cavity-QED system, 
where |+) denotes the excited QD state, |— ) denotes the 
ground state, |1) and |0) represent the first excited and ground 
state of the cavity mode. Also shown is the phonon reservoir 
as blue (multiple) lines. The terms {B)g and {B)rjx represent 
a phonon- modified coherent reduction {{B) (T) < 1) in the 
exciton- cavity coupling rate and the exciton pump rate, re- 
spectively (see text for details), (c) Schematic of a dot-driven 
cavity-QED system with same energy level notation as in (b). 



pump frequency, cjl, we first introduce the model Hamil- 
tonian describing a cavity-QED system where the QD 
interacts with an acoustic phonon reservoir: 



where 6g(6j) are the annihilation and creation operators 
of the phonon reservoir, a is the leaky cavity mode an- 
nihilation operator, (j+ (annihilation) and cf~ (creation) 
are the Pauli operators of the electron-hole pair or ex- 
citon; Ac^L = uoa — {oi = x^c) are the detunings of 
the exciton {uOx) and cavity {uOc) from the coherent pump 
laser (a;^), and g is the cavity-exciton coupling strength. 

The pump term, i^^/^^e' accounts for the coherent drive 
on the cavity-QED system; for a QD (exciton) driven 
system, i^Jrive — ^Vx{o'~^ -^^~)^ while for a cavity driven 
system, i^^^ive ~ ^^c(^ + ^^). The defined pump rate, 
r]x/ci is two times the classical Rabi frequency. 

Transforming to the polaron frame, we eliminate the 
QD-phonon coupling and introduce a renormalized dot- 
cavity coupling strength^^. For the case of the QD 
driven system, the polaron transformation also results 
in a renormalized Rabi frequency, defined below. The 



polaron transformation^^ can be written as 
= exp(5')i^exp(-5'), 

where 
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The transformed Hamiltonian becomes 

H'^y^ = h{AxL - Ap)(7+(7- + hAcLa^a + {B)Xg, (4a) 
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where flq = (o^qhq) = [e^^^-? — is the mean phonon oc- 
cupation number (Bose-Einstein distribution) at a bath 
temperature, T = l/k^p. For clarity, we will henceforth 
assume that the polaron shift is implicitly included in our 
definition of ujx (one should, however, keep in mind that 
this shift is temperature dependent). For a dot (exciton) 
driven system, Xg and Xu are defined through 

Xg = hg{a^a~ + a~^a) + hrix{cF~ + a^), (8a) 
Xu = ih[g{cF^a — cr~) + ihirix{cr^ — a~)], (8b) 
and for a cavity driven system, 

Xg = hg{a^a~ + a^a), (9a) 
Xu = ihg{d-~^a — d^a~). (9b) 

Worth to note is the slightly unusual definition of the 
system Hamiltonian, Eq. (|4a|) . The usual (but in gen- 
eral, incorrect) decomposition of the system Hamiltonian 
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to include only the noninteracting QD and cavity parts 
does not take into account the effect of the coherent cw 
drive on the system Hamiltonian. As the cavity and 
the QD systems are internally coupled, as discussed by 
Carmichael and Walls^^, this leads to violation of detailed 
balance. The system Hamiltonian written above leads to 
the correct form of the density operator while preserving 
detailed balance. Moreover, it includes the effect of dot- 
cavity coupling and the dot-cw driving on the coherent 
part of the Hamiltonian to all orders. 

Next, we unitarily transform to a frame of reference 
defined by this system Hamiltonian which we will use 
to obtain a time-convolutionless ME; the net effect of 
this transform, e.g., in the case of resonance fluorescence 
through an exciton driven system, results in Mollow 
triplet peaks that sample the asymmetric phonon bath at 
the dressed eigenfrequencies — as determined by the mod- 
ified system Hamiltonian. Somewhat similar techniques 
(i.e., bath sampling at the dressed resonances) have been 
employed to study atomic dynamics in generalized {en- 
gineered) photon reservoirs, including photonic band gap 
materials^ and squeezed reservoirs^^. 

Phenomenologically, we include the radiative decay of 
the QD and the cavity mode as Liouvillian superopera- 
tors acting on the reduced system density matri:^. In 
addition, we incorporate an additional pure dephasing 
process beyond the IBM with a rate 7' — this accounts 
for the broadening of the zero-phonon line (ZPL) with 
increasing temperature o^^i^^ — . Though there is some 
controversy about what causes the broadening of the 
ZPL, e.g., spectral diffusion, anharmonicity effects^i^ 
phonon scattering from interfaces and a modified 
phonon spectrum^^ , it is well known that the ZPL broad- 
ens as a function of temperature with a Lorentzian scat- 
tering process; thus, we treat the broadening of the ZPL 
phenomenologically, while accounting for broadening as 
a function of temperature similar to experiment a^^i^^i^^ , 
with y{T) scaling as ~ IjiieV/K. The various superop- 
erators act on the reduced system density matrix , and 
are defined through 

L{p) = —{2a~pa~^ — cr~^a~p — pcr~^a~) 
+ t<i(2apa^ — a^ap — pa^a) 

i 

+ — (2(Jiip(Jii - aiiaiip - panaii) , (10) 

where 2i<i is the cavity decay rate, 7 = is the ra- 

diative decay rate, and an = cf^cf~ . The radiative de- 
cay rate has an additional renormalization by a factor of 
(5)^, which reduces the effective radiative decay rate in 
the presence of phonons^ . 

We then derive a time convolutionless ME for the re- 
duced density operator, p(t), of the cavity-QED system"^^ 
in the second-order Born approximation (for incoherent 
bath coupling). The time-convolutionless form of the 
ME, though local in time, is known to capture non- 
Markov effects due to the reservoir—. However, for our 
analysis, we will make a Markov approximation as typical 



phonon processes are substantially faster (i.e., a few ps) 
than the relevant system dynamics by at least an order of 
magnitude. This allows us to obtain effective rates which 
naturally depend on the spectral densities of the phonon 
spectral function that are locally sampled by the dressed 
resonances. We have checked that the Markov limit of 
the fully non-Markovian time-convolutionless ME is rig- 
orously valid for the system and excitation (cw) cases of 
interest. Thus, while it is straightforward to carry out 
non-Markov calculations, it is not necessary here — they 
give identical results. 

In the interaction picture described by H'^y^^ we con- 
sider the exciton-photon-phonon coupling H[^^ to second 
order (Born approximation), and trace over the phonon 
degrees of freedom to obtain a Markovian time convolu- 
tionless ME^^^^^: 

m—g,u ^ 

Xra, e-*^->^/^X™e^^=-^/V(t)] + H.c)i , (11) 

where Gg/u{t) = ((<?/« (OCg/«(0))- The polaron Green 
functions are^i^, 

Gg{t) = {B)^{cosh[(t>{t)]-l), (12a) 
G„(i) = (B)2sinh[0(i)], (12b) 

which depend on the phonon correlation function, 

roo J((jj) 

(j){t) = / duo^Y- [coth(/3ncj/2) cos(a;t) - i sin(cjt)] , 
Jo ^ 

(13a) 
(13b) 



where J{w) is the the characteristic phonon spectral 
function, defined in this work as 



J{uj) = ap uj'^ exp I — 



(14) 



This form of the spectral function [Eq. (fT4|) ] describes 
the electron-LA(longitudinal acoustic)-phonon interac- 
tion via a deformation potential coupling which is the 
main source of dephasing in self- assembled InAs/GaAs 
QDs. For all our calculations that follow, we use param- 
eters suitable for InAs/GaAs QDs^^, with (^5 = 1 meV 
{(jJi) is a high frequency cutoff proportional to the inverse 
of the typical electronic localization length in the QD) 
and ap/(27r)^ = 0.06 ps^; these values vary somewhat in 
the literature, though we have taken ours from fitting 



recent experimentsi 



24.39.40 



Using the parameters above. 



e.g., at T = 10 K, yields a polaron shift, Ap = 42 peV ^ 
and a Franck-Condon renormalization, {B) = 0.84. With 
these phonon parameters, we already see that clearly the 
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coherent renormalization effects will be important for an- 
lyzing PL intensity for QD-cavity systems, even at rela- 
tively low phonon bath temperatures. 

We briefly mention that there are other electron- 
acoustic phonon scattering models that can go beyond 
the polaron ME approach. For example, McCutcheon et. 
als^ recently introduced a more general ME technique 
to describe the non-equilibrium dynamics of a QD sys- 
tem interacting with a phonon reservoir based on a vari- 
ational formulation (with no cavity coupling). This ele- 
gant approach extends the validity of the ME to parame- 
ter regimes, rjx ^ ^6, where the ME in the polaron frame 
can break down. However, the pump parameter regimes 
that we study in this work {r]x/c <^ ^b) are well within the 
domain of validity of our ME, so we can safely use the po- 
laron ME, while also accounting for cavity coupling^^. A 
more general description of the system dynamics valid in 
all regimes can be obtained using a quasi-adiabatic path 
integral approach^^^^^. The beneflts of our polaron ME 
is that the solution, even with multiphonon and multi- 
photon effects included, is relatively straightforward, and 
it has already been used to help explain experiments for 
coherently-excited dots in the regime of cavity- QEE^i^. 
Moreover, as we will show below, one can derive a user- 
friendly Lindblad ME that contains many of the key fea- 
tures of phonon interactions in cavity- QED systems. 

For numerical calculations, we solve the above ME with 
steady-state pumping (i.e., r]x/c are time-independent), 
with the exciton initially in the ground state. Prior to 
these dynamical calculations, we compute the phonon 
scattering terms in Eq. (pTj) . whose solution is natu- 
rally problem-dependent (through H'^y^). Thus there 
are no flxed phonon scattering rates for analyzing QD 
power broadening as a function of pump power, as the 
phonon scattering rates are pump-dependent. The same 
arguments apply for studying power broadening as a 
function of temperature; one must obtain the phonon- 
induced scattering rates for each pump value and tem- 
perature. Experimentally, the intensity PL lineshape 
is usually obtained by measuring the QD exciton in- 
tensity {Ix) or cavity mode intensity (/c) as a function 
of increasing pump fleld. To connect to these quanti- 
ties, we solve the above ME in a Jaynes-Cummings basis 
with states |0) , , \IU) , \2L) , \2U) , • • • , and compute 
the steady-state exciton and cavity photon populations, 
fix = (cr+cr~) \ss oc Ix and fic = {a^a)^^ oc /c. Deflning 
the photon states |n), with n = 0,1,2,---, and exciton 
states 1+/—), then the Jaynes-Cummings ladder states 
are related to the bare states, e.g., through: |0) = |— ) |0), 
lli) = ^(h) |1)- 1+) |0)), \1U) = ^(1-) |1) + |+) |0)). 

In our calculations we make partial use of the quan- 
tum optics toolbox by Tan^^, and flnd that truncation 
to two-photon-correlations (two photons or 5 states) is 
sufficient /necessary for all the dot driven simulations, 
while truncation to six-photon-correlations (six photons 
or 13 states) is sufficient /necessary for all cavity calcula- 
tions that follow. The role of multiphoton effects de- 
pends on the value of the dot-cavity coupling rate 



which we choose to be ^ = 20 /ieV — consistent with typ- 
ical semiconductor cavi ty-QE D power broadening exper- 
iments (e.g., see Refs. |3Q||35[ ). A detailed discussion of 
the role of multiphoton, and multiphonon processes, is 
presented in Sec. IIVBI 

III. EFFECTIVE PHONON MASTER 
EQUATION OF THE LINDBLAD FORM 

Our polaron ME [Eq. (pT]) ] includes both coherent and 
incoherent contributions from electron-phonon scatter- 
ing, but some care and insight is needed in extracting 
the relevant incoherent scattering rates. It is therefore 
instructive to construct a simplifled phonon-modifled ME 
of the Lindblad form, which we call an effective phonon 
master equation (EPME); we do this by simplifying the 
term, e~'^^^y^^^^Xme'^^^y^^^^, appearing in the fuh time- 
convolutionless ME [Eq. (pT]) ]. The resulting Lindblad- 
form ME enables a very simple numerical solution and 
facilitates the extraction of various phonon-induced scat- 
tering rates in a clear and transparent way. We expect 
that the integral in Eq. (fTTj) can be approximated, un- 
der certain circumstances, by only including the phase 
evolution of the operators Xg^u with respect to the non- 
interacting part of the system evolution. Further, for a 
QD-driven system we only include terms proportional to 
g'^ and r]^ and ignore cross terms proportional to grjx'', the 
inclusion of the cross terms do not preserve the Lindblad 
form and contribute very little to the overall broaden- 
ing lineshape as can be demonstrated numerically. For a 
cavity-driven system, we again include the phase evolu- 
tion of the operators Xg^u with respect to the noninter- 
acting part of the system evolution; however, the effective 
Lindblad description has only contributions which are 
proportional to g'^ — since Xg^u do not depend on rjc We 
will, of course, compare the EPME solution with the full 
numerical solution of the polaron time-convolutionless 
ME, i.e., Eq. ([11]); the prime purpose of the EPME is 
to help elucidate the physics of phonon-induced incoher- 
ent scattering, though we will highlight regimes where it 
can work quite well in accurately describing the full char- 
acteristics of the entire power-broadened PL hneshape. 

We postulate that the dynamics of the QD driven sys- 
tem can now be approximately described through 

g = ^[H:^,,p{t)] + L{p) + LM (15) 
where Lph{p) {^ph^ refers to phonon) is given by 

+ ^L{a+a) + ^L{a^a-), (16) 

and the superoperator L{D) is deflned as 

L{D) = 2Z)pi3t - DWp - pDW. (17) 
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The above effective ME [Eq. (p!5]) ] has a remarkably sim- 
ple form, and its general format should be familiar to 
many researchers who have been using atomic cavity- 
QED models to connect to experimental data using semi- 
conductor cavity- QED systems. However, it must be 
used with caution, as it is only valid within certain 
regimes where the above noted approximations are good. 
The phonon- mediated rates, which drive the effective 
Lindblad dynamics, are derived to be (see Appendix [A]) : 



dre 



(18) 
) . 

(19) 



where /S^cx = O0c — (^x is the cavity-exciton detuning. Fig- 
ure [2] shows a schematic of the various effective phonon- 
scattering processes: F^^ describes phonon-assisted inco- 
herent excitation and EID (pump-induced broadening); 
F^^ described enhanced radiative decay and EID; F^^^ 



describes (the somewhat unlikely scenario of) exciton ex- 
citation via the emission of a cavity photon, and F^^^ 
describes the process of cavity excitation (cavity feeding) 
via the absorption of a photon. Importantly, all of these 
scattering events are driven by electron-phonon interac- 
tions and they cause effects that are significantly different 
to simple pure dephasing models. In fact, pure dephas- 
ing through 7' only results in Lorentzian coupling, and 
is found to play a minor role in what follows below. 
Our formalism above shows that for a cavity-driven 

system, F^^ =0 and A^^ =0, and there is no 
phonon-induced EID due to the lack of any coupling of 
the drive with the phonon reservoir. We therefore expect 
(and find) 

substantially different intensity power broadening be- 
tween dot-driven and cavity-driven systems; both of 
these exciton-driven and cavity-driven models are also 
markedly different to simple atomic models. We high- 
light that a similar exciton-cavity (feeding) rate, F^J^^ , 
has been derived by Xue et al?^ and by Hohenester— , 
though these were obtained for an undriven cavity- QED 
system. Both of these useful approaches also use a po- 
laron frame to describe the incoherent scattering, though 
the end equations have some potential problems for small 
cavity-exciton detunings (where, admittedly, these effec- 
tive rates are at best approximate anyway), and neither 
approach includes a coherent (temperature-dependent) 
reduction in ^ ^ {B)g. For example, Hohenester— de- 
rives the following exciton-cavity feeding rate: F^^^ = 



2^^ Re [Jq^ (ir e~*^^^'^e~'^^'^^] , which has a similar form 
to our Eq. ([19]) (apart from the sign of the phase and 
the need to subtract of a background term for small de- 
tuning), but does not include (5(T))^ — an important 
temperature-dependent term. The general predictions 
using this rate formula are consistent with experiments^. 
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FIG. 2: (Color online) Schematic of the phonon-scattering 

processes: (a) T^t , (b) V;' , (c) r^>, and (d) T^I"", for 
a QD-driven system introduced in the effective Lindblad de- 
scription and defined in Eqs. (|18p - (|19p . The various symbols 
are the same as in Fig. [1] Processes (a) and (b) describe 
phonon-induced incoherent excitation and excitation induced 
dephasing; process (c) describes phonon-induced exciton exci- 
tation via cavity decay, though this process has a small prob- 
ability because of relatively large cavity decay rate (n); pro- 
cess (d) describes phonon-induced cavity excitation through 
the decay of an exciton; at low temperature this (cavity feed- 
ing) process occurs primarily via phonon emission, and will 
be more efficient for a cavity mode that is red shifted with 
respect to the exciton. The sign of the cavity detuning in (c- 
d) is just an example, and note that the phonon interactions 
with A cx — ^ or l\cx — 

A will be different (especially at 

lower temperatures). 



and all these aforementioned polaron formalisms pro- 
duce qualitatively the same trend as a function of cavity- 
exciton detuning (compare results in Refs. [25II27I with 
those in Figs. EE]). 

An alternative effective Lindblad ME with the same 
process identified in Refs. [25II27I and above (i.e., for 
exciton-cavity coupling), was recently presented by Ma- 
jumdar et. ali^ and used in part to study the role of 
phonon scattering for the cavity-emitted resonance fiu- 
orescence spectrum. The infiuence of phonons was in- 
cluded as two additional incoherent decay terms which 
were to second-order in the QD-phonon coupling. Unlike 
the polaronic approaches above, electron-phonon inter- 
actions were included only to first order, which is gen- 
erally not valid in these cavity- QED systems — even at 
low temperature^^ (see Sec. lIVBp . More problematic 
is the fact that the effects of the coherent drive on the 
phonon reservoir and the associated EID effects are miss- 
ing; in contrast, we find these to be the dominant source 
of broadening from electron-phonon scattering. The need 
for EID processes in coherently driven semiconductor 
cavity- QED system has already been shown for the Mol- 
low triplet, both experimentall}^ and theoretically^. 

In addition to the phonon-induced Lindblad decay 
rates above, one also has phonon- mediated frequency 
shifts beyond the polaron shift. The effective Hamilto- 
nian, describing the coherent part of the system evolu- 
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tion, -ff®?,, becomes 



(20) 



with 



dr 



Acr"'"a/a^cr / 7->\ 2 2 t 



(21) 



(22) 



where A^^^, A^J^^ , A^^ and A^^ are the Stark shifts 
(which scale proportionahy with {B) g'^ or {B) r]^). 



IV. NUMERICAL RESULTS 

A. Role of Phonon Scattering on Intensity Power 
Broadening: Effective Phonon ME Versus the Full 
Time-Convolutionless ME 

We first investigate the role of the four phonon Lind- 
blad terms in a typical power-broadened intensity PL 
computed with our EPME [Eq. ([15])] and compare with 
the full solution (Eq. (pT]) : time-convolutionless ME). 
The main parameters are listed in Fig. [3l and we have 
adopted system parameters and coupling constants sim- 
ilar to those in recent semiconductor experiments^Ui^. 

In Fig. [3] we study the role of F^^, F^^, F^J^^" and 
F^^^, as defined in Eqs. (jl8til9p . on the power broaden- 
ing lineshape {Ic oc nc); here we use rjx = 40 /ieV for a 
QD driven system at a bath temperature of T = 4 K, 
and study two different cavity-exciton detunings, (a-b) 
Acx = 3 meV and (c-d) Acx = —3 meV. The corre- 
sponding peak Uc that results from this interaction is 
around 4 x 10^. To better highlight the various scattering 
mechanisms, we include only one of the Lindblad terms 
in each calculation, as labeled in the plots. By looking 
at Fig. [3Ka,c), it is clear that the process L{d-~^) is pri- 
marily responsible for incoherently exciting the phonon 
sidebands [cf. [3](a)] and EID; while process L{a~) in- 
troduces further pump-dependent EID, as will be high- 
lighted in detail later (note that this process has the 
same Lindblad operator terms as 7). The broad back- 
ground centered at cjl — = 1 meV is present only for 
the case of F^ . The exciton-cavity scattering processes, 
L{d^d-~) and L((j+a), account for cavity excitation and 
cavity destruction, respectively, by phonon-assisted pro- 
cesses and these affect the relative magnitudes of the 
cavity measured intensity PL at different temperatures 
and drives. Figures [3jb,d) demonstrate that L{d^a~) 
is the main cavity-exciton coupling (feeding) term; this 
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FIG. 3: (Color online) (a) Influence of the effective phonon 
scattering rates, Tp^ (red, lower solid line) and F^^ (black, 
dashed line) [defined in Eq. (|18|) ]. on the intensity PL, /c, 
for a QD driven system at T = 4 K with Acx = 3 meV. 
Also shown is the full polaron ME solution (grey, upper solid 
line) [Eq. (|lip ]. We plot the cavity intensity (Ic) for exciton 
excitation as a function of QD-laser detuning and show the 
contribution of one effective Lindblad rate per calculation; 
the collective infiuence of these processes is shown later in 
Fig. [5] and compared with the full solution, (b) As in (a), but 

for effective phonon scattering rates, Tpl^ (red, lower solid 

fine) and F^^^ (black, dashed line) [defined in Eq. The 
system and material parameters are: j — 2 /xeV, k — 50 /xeV, 
g = 20 /xeV, r]x = 40 /xeV, 7^(2 K) = 2 /xeV, and we compute 
(B) (4K) = 0.91. (c-d) As in (a-b) but with Acx = -3 meV. 
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FIG. 4: (Color online) (a) Infiuence of the role of the effective 

phonon scattering rates, F^^^ (red, lower, dashed line) and 

Tp^^ (blue, upper, dashed line), on the intensity PL, /x, for 
a cavity driven system at T = 4 K for Acx = —0.5 meV. 
Also shown is the full polaron ME solution (solid grey lines) 

[Eq. (pT]) ]. For clarity, the Tp^^ curves are shifted vertically 
by 0.3, along with the full polaron ME solution. We plot the 
contribution of one effective Lindblad rate per calculation, 
(b) As in (a) but for Acx = 0.5 meV. The parameters are the 
same as in Fig. O but with ric = 40 /xeV (cavity pumping). 



mechanism results in enhanced cavity photon numbers at 
the exciton transition, especially when the cavity is red 
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FIG. 5: (Color online) Phonon bath at T = 4 K. (a-b) Nor- 
malized cavity mode intensity (/c) for a dot-driven system as 
a function of QD-laser detuning for two different dot-cavity 
detunings, (a) Acx = 3 meV and (b) Acx = —3 meV, and for 
two different values of the cw laser Rabi frequency (orange, 
lower solid line, corresponds to r]x = 20 /xeV; and grey, upper 
solid curve, corresponds to rjx = 40 /xeV). Also shown (black, 
dashed lines) are the intensity PL obtained using the effective 
Lindblad form of the full time-convolutionless ME. Note that 
we have vertically shifted the intensity PL for r/x — 40 /xeV 

by 0.2 for clarity, (c) Plot of the phonon rates F^^ (see 
text) for Acx = —3 meV as a function of QD-laser detuning. 
Note (-Acx) = F'^ti^cx). (d) Plot of the phonon rates, 
^ph ^ , as a function of QD-cavity detuning. 



shifted from the exciton — since phonon emission is favor- 
able at lower temperatures. In contrast, the L((j+a) pro- 
cess gives no noticeable exciton-cavity coupling because 
of the fast cavity decay rate^i. 

In Fig. |4] we carry out a similar exercise for a cavity- 
excited system {r]c = 40 /ieV), calculating Ix (oc nx)^ 
where we study the influence of L{a^a~) and L{a~^d) on 
power broadening; here we find excellent agreement with 
only the F^J^^ scattering term (lower red, dashed curve) 
compared to the full ME solution (lower grey, solid line), 
which results in a significant exciton-cavity feeding pro- 
cess via phonon emission [cavity is red detuned in (a), 
cf. Fig. Eld)]. Again we find that L{a~^a) gives no no- 
ticeable cavity feeding. For this cavity driven system, we 
have chosen Acx = ±0.5 meV instead of Acx = ±3 meV 
(which we chose earlier for the exciton driven system); 
this is because a strong exciton-driven system invariably 
kicks up the phonon sidebands even at low temperatures 
that can swamp the emission at the cavity mode; so we 
use a larger cavity-exciton detuning for the exciton driven 
case. The cavity driven system is thus much cleaner 
to analyze for smaller cavity-exciton detunings, and the 
exciton-measured intensity PL is also substantially re- 
duced for larger detunings. Note that the corresponding 
peak fix for this excitation regime is around 0.3. 



FIG. 6: (Color online) As in Fig. [5]but with the phonon bath 
at T = 20 K. The system parameters are identical to those 
given in Fig. [3l except that 7^(20 K) = 20 /xeV and we now 
compute {B) (20 K) = 0.73 (cf. (B) (4K) = 0.91). 



In Figs. [5Ka,b) we plot the normalized cavity mode 
intensity, Ic cx nc, for a dot-driven system as a func- 
tion of QD-laser detuning, again for the two different 
cavity-exciton detunings, (a) Acx = 3 meV and (b) 
^cx = — 3 meV, at T = 4 K. We also show the inten- 
sity PL for rjx = 20 /ieV and rjx = 40 /ieV. To study 
the effects of increasing temperatures, in Fig. [6] we plot 
the /c at T = 20 K. In these graphs, we show the total 
power broadened intensity PL obtained using the EPME 
[dashed lines in (a-b)] and compare with the full polaron 
ME solution [solid lines in (a-b)]; given the approxima- 
tions made in the derivation of the EPME, the agree- 
ment is remarkably good. The frequency-shift terms, 
A^^, are found to be very small here and can be ne- 
glected for the cases shown. Note that the phenomeno- 
logical pure-dephasing rates are temperature dependent, 
where we choose y(4K) = 2 /ieV and 7^(20 K) = 20 /ieV 
(e.g., see Refs. |2Q||53[ ). The thermal expectations of 
the phonon displacement operators are calculated to be 
(5)(4K) = 0.91 and (5)(20K) = 0.73. These re- 
sults suggest that the dynamics can be well described 
by our EPME, by essentially only including three sepa- 
rate phonon scattering effects — since, from the findings 
above, L((j+a) can be safely neglected. 

In Fig. [51(c) we plot the corresponding phonon scatter- 
ing rates, F^^ and F^^ , as a function of QD-laser detun- 
ing for r]x = 40 /ieV, at T = 4 K; in Fig.[6^c) we plot these 
scattering rates for T = 20 K. Since the rates depend on 
QD-laser detuning and the pump strength, they are obvi- 
ously important for understanding power broadening in 
a cavity-QED system. In Figs. [5]-6(d) we plot the rates, 
F^^^ and F^^^ , as a function of QD-cavity detuning for 
two different temperatures; the regions marked by the 
vertical lines indicate the chosen detunings in Figs, (a) 
and (b) — so note that these particular rates are fixed as a 
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FIG. 7: (Color online) (a) and (b) Relative QD intensity for 
a cavity- driven system as a function of QD-laser detuning for 
Accc = — 0.5 meV and phonon bath temperature (a) T = 4 K 
and (b) T = 20 K and for two different values of the cw 
laser Rabi frequency (orange, lower, solid line corresponds to 
r]c = 40 /ieV, and grey, lower, solid line corresponds to r]c = 
60 /ieV). Also shown (black dashed curves) are the intensity 
PL obtained using the effective Lindblad form of the full time- 
convolutionless ME. The parameters were as follows: 7^ = 
2 //eV, = 50 /ieV, 7' = 2 /xeV at T = 4 K, and 7' = 20 /xeV 
at T = 20 K. (c) Same as in (a) and (b) but with Acx = 
0.5 meV. 



function of QD-laser detuning. As discussed earlier, the 
scattering term, describes a process which involves 

de-exciting the QD and exciting the cavity mode, aided 
through phonon emission or absorption. 

The collective influence from the various phonon scat- 
tering terms, discussed above, results in broadening 
(EID) of the QD exciton resonance, incoherent excita- 
tion of the phonon bath, and significant exciton-cavity 
coupling (or feeding); the first two of these are pump- 
dependent, through ~ (5)^7^^, while the latter process 
scales with (B)'^ q'^. The trends of the exciton-cavity 



{B)y_. 

feeding rate, T'^^^ , are consistent with the results of 
Refs. I25lf26. where one observes a peak scattering rate 
at around A 

cx J- — 

2 meV (depending upon the tem- 
perature). Note that at T = 4 K and Acx = — 3 meV, 
per a dominates, whereas F" ^ is much larger for Acx = 
3 meV. However, at T = 20 K, the two rates are much 
closer to each other for the different detunings as the de- 
pendence of the rates on QD-cavity detuning becomes 
more symmetric with increasing temperature^^. 

It is interesting to note that the broadening of the 
QD exciton resonance as a function of QD-laser detuning 
closely mirrors the PL lineshape associated with the lin- 
ear exciton spectrum obtained using the IBM. A similar 
observation was demonstrated by Ahn et. al.^^ ^ where 
the effects of electron-phonon coupling in QDs (with no 
cavity) on nonstationary resonance fluorescence spectra 



were studied; the resonance fluorescence dynamics of the 
QD electronic transition was shown to have a strong 
dependence on the duration of the laser field, and by 
increasing duration of the laser pulse, the background 
phonon continuum was strongly excited. 

Similar to the QD driven system, in Fig. [3 we also 
investigate the exciton PL characteristics, Ix oc fix, but 
with cavity excitation (investigated in more detail in Sub- 
sect. |IVC]). As anticipated from Fig.lH we again obtain a 
very good fit between our EPME results and the full po- 
laron time-convolutionless ME solution. In general, the 
effective phonon Lindblad solution is expected to closely 
mimic the full solution here for large QD-cavity detun- 
ings (|Aca;| ^ g^Vc)^ though we find that it can work well 
over a wide range of excitation conditions. 

A preliminary analysis of the role of the various scat- 
tering processes on the intensity PL is presented in Ta- 
bles Hl lVIl Here we focus on the qualitative influence of 
the various scattering rates on the exciton (x) and cav- 
ity (c) FHWM broadenings, but later we will investigate 



TABLE I: Dependence on the x/c FWHM of the intensity PL 
on 7, for a QD-driven cavity measured PL (x FWHM, with 
Acx = 3 meV), and for a cavity-driven QD measured PL 
(c FWHM, with Acx = 0.5 meV). The QD/cavity parameters 
are z^: = 50 /xeV, g = 20 /xeV, //(cc/c) = 30 /xeV, and 7^ = 2 /xeV 
(4 K). No coupling to phonons is considered and all numbers 
are in units of /xeV. 







X FWHM (QD driven) 


c FWHM (cavity driven) 


7 


= 2 


120 


107 


7 


= 4 


104 


107 


7 


= 6 


98 


107 



TABLE II: Dependence on the x/c FWHM of the power 
broadened lineshape on for a QD-driven cavity measured 
system (x FWHM, with Acx — 3 meV), and a for cavity- 
driven QD measured PL (c FWHM, with Acx = 0.5 meV). 
The QD/cavity parameters are the same as in Table HI and 
all units are in /xeV. 





X FWHM (QD driven) 


c FWHM (cavity driven) 


K = 10 


120 


41 


= 30 


120 


60 


= 50 


120 


107 



TABLE HI: Dependence on the x/c FWHM of the power 
broadened lineshape on 7^, for a QD-driven cavity measured 
PL system (x FWHM, with Acx = 3 meV), and for a cavity- 
driven QD measured PL (c FWHM, with Acx = 0.5 meV). 
The QD/cavity parameters are the same as in Table m and 
all units are in /xeV. 







X FWHM (QD driven) 


c FWHM (cavity driven) 


7' 


= 2 


120 


107 


7' 


= 4 


146 


107 


/ 

7 


= 6 


169 


107 
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TABLE IV: Dependence on the x FWHM of the power broad- 
ened intensity PL of a QD-driven cavity measured system 
(Accc = 3 meV), on the two phonon Lindblad processes, F^/^ 

and Fp/j , at T = 4 K, for various values of r\x. The general 
system parameters are the same as in Table [U Also consid- 
ered is the case of no incoherent phonon coupling with {B) — 1 
(i.e., also no coherent phonon coupling), and (B) = 0.91 at 
T = 4 K (with coherent phonon coupling). The inclusion 
of l\B) at 4 K here is to make a better comparison with the 
same effective r\x (^ (5) r\x) and effective g (^ (B) ^); in re- 
ality, the no-phonon case will be different because of the lack 
of coherent phonon renormalizations (later we will show this 
difference more explicitly, e.g., in Figs. [TQ1IT3)) . All numbers 
are in units of /ieV, and note that the broadening here is in- 
dependent of Fp^ " and . The last column gives the full 
polaron ME result [Eq. (ITTl) ] 



No 

(B) = l 



ry^ 20 80 
ry^ = 40 159 
T]^ = 60 240 



No F^;, 


J- vh 


per 


EPME 


Full 


(5) (4 K) 






73 


71 


72 


70 


72 


146 


138 


139 


134 


151 


219 


251 


200 


223 


390 



TABLE V: Dependence on the c FWHM of the power broad- 
ened intensity PL of a cavity-driven QD measured system. 



and r^l" 



at 



on the two phonon Lindblad processes, F^;^ ' 
T = 4 K, for various values of g. Also considered is the case 
of no phonon coupling, with two {B) values as above. The 
general system parameters are the same as in Table H and all 
numbers are in units of /ieV. 



20 
40 
60 



No 

(S) = l 

101 
101 
101 



No F-^ 
(5) (4 K) 

101 
101 
101 



ph 



101 
101 
101 



ra^cr 
ph 



102 
109 
115 



TABLE VL Dependence on the c FWHM of the power broad- 
ened intensity PL of a cavity-driven QD measured system 
(Acx = 0.5 meV), on the two Lindblad processes, Tp^^^ and 



- ph 



at T = 4 K, for various values of r]c. Also considered 
is the case of no phonon coupling with two (B) values. The 
general system parameters are the same as in Table H] and all 
numbers are in units of /xeV. 



the full power-broadened PL curves in detail. In Table |T] 
we study the dependence of the x/c FWHM of the in- 
tensity PL on 7 (the radiative decay rate). Increasing 7 
results in reduced broadening of a cavity-measured QD 
driven PL lineshape (x FWHM). Increasing 7 broadens 





No 


No r^^ 


J- ph 


-L ph 


EPME 


Full 




{B) = l 


{B){4 K) 










= 20 


100 


100 


100 


100 


100 


100 


= 40 


101 


101 


101 


103 


103 


103 


= 60 


102 


102 


102 


107 


108 


108 



the QD resonance resulting in reduced coupling between 
the cavity mode and the QD and subsequent narrowing 
of the cavity- measured QD driven PL lineshape. How- 
ever, 7 has little effect on the power broadening of a 
QD-measured cavity driven intensity PL (c FWHM). In 
Table im we present numerical results for the x/c FWHM 
of the intensity PL for different values of k,. Increasing k 
here has no observable effect on QD-driven cavity mea- 
sured intensity PL, but increasing tz results in increased 
broadening of a cavity-driven QD measured intensity PL. 
In Table Hill we study the dependence on 7^, i.e., pure de- 
phasing of the ZPL. Increasing 7' increases broadening 
of QD-driven cavity measured intensity PL, but has neg- 
ligible effect on the cavity driven QD-measured system. 

In Table HVl we study the role of the two phonon Lind- 
blad processes, F^^ and F^^ , on a QD- driven cavity- 



measured system at T = 4 K, for various values of rjx; 
the other two phonon processes, i.e., T^^^ and F^J^^ , 
do not influence the power broadening lineshape here. 
As can be seen, increasing the influence of these effective 
Lindblad processes by increasing r]x reduces the linewidth 
of the cavity measured intensity PL. We also see that the 
X FWHM broadening for increasing drives may be dras- 
tically underestimated by the effective phonon ME, even 
though the qualitative trends of the full intensity PL are 
similar; this is because we have neglected the influence 
of the drive on the phase integrations that enter the full 
incoherent phonon scattering [see Eqs. ([8a|) and (pT]) ]. 

In Table El we list the c FWHM of the power- 
broadening lineshape of a cavity-driven QD measured 
system on the two Lindblad processes, T'^^^ and F^J^^ , 
at T = 4 K, for various values of g. The Lindblad process 
caused by F^^ ^ has no effect on the c FWHM. However, 

increasing F^J^^ (the cavity feeding process) results in 
some increased broadening. Finally, in Table IVll we show 

the influence of these phonon processes {T^j^ ^) with 
increasing ?7c, which also show a slight increase of the c 
FWHM with increasing drives. Now we see that the c 
FWHM broadening for increasing drives is extremely well 
reproduced with the effective phonon ME; this is because 
the cavity pump no longer enters H^^y^ [see Eqs. ([9a|) and 
(pT]) ]. In both Tables [VTlVIi we see no role from the slight 
coherent reduction of {B)g at 4 K. In general we see that 
driving via the exciton or the cavity can yield drastically 
different results. 

To summarize this subsection, we have identified three 
main Lindblad phonon processes that contribute to the 
intensity PL of a cavity-QED system. We have also found 
earlier that the effective scattering rates associated with 
these processes can be calculated, under certain detuning 
conditions, from simple analytical solutions [Eqs. [T8lfT9] : 
these effective Lindblad solutions are compared to the 
full polaron ME solution and found, in certain cases, to 
yield very good agreement. In this way we can also argue 
the underlying physics of the identified phonon scattering 
processes. However, there can be noticeable differences, 
especially for the measured x FWHM broadening of an 
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FIG. 8: (Color online) (a-b) Cavity mode intensity PL for a 
QD driven system {rjx = 30 /xeV) as a function of QD-laser 
detuning for a one-phonon theory (red dashed line) and the 
full polaron model (grey solid line) with Acx = —3 meV, 
for phonon bath tempereatures (a) T = 4 K and (b) T — 
20 K. (c-d) Relative QD (exciton) intensity for a cavity-driven 
system {rjc = 30 /xeV) as a function of QD-laser detuning for 
a one-phonon theory (red dashed line) and the full polaron 
model (grey solid line) with Acx = —0.5 meV, for phonon 
bath temperatures (c) T = 4 K and (d) T = 20 K. 



exciton driven system. 

In what follows below, we will use the full polaron ME, 
and first verify the general need for for multiphonon and 
mutiphoton effects. 



B. Influence of Multiphonons and Multiphotons on 
the Intensity PL 

Our time-convolutionless ME above [Eq. (fTT]) ] utilizes 
the polaron frame which allows for a nonperturbative 
treatment of phonons. This enables one to use the full 
IBM machinery to compute the phonon correlation func- 
tions. It is also useful to look at the one phonon limit of 
the time convolutionless ME by expanding the phonon 
Green functions to lowest order in the phonon coupling. 
In this limit we can expand the phonon correlation func- 
tion 0(t) to lowest order in the dot-phonon coupling con- 
stant as follows: Gg{t) and Guif) (j){t) where we 
have used {B) 1. This then connects to a weak-phonon 
coupHng (i.e., perturbative) approach. 

In Figs. IHl^a) and (b) we plot the cavity mode inten- 
sity (/c) for an exciton driven system as a function of 
exciton-laser detuning for a one phonon and the full po- 
laron (multi-phonon) solution. We recognize that the 
one phonon scattering process tends to deviate from the 
full polaron intensity PL, especially noticeable at higher 
temperatures and for larger driving field strengths. Thus 
even for low temperatures, a weak-phonon coupling the- 
ory can break down. We further find that the one phonon 



FIG. 9: (Color online) (a-b) Exciton intensity PL for a cavity 
driven system {rjc = 30 /ieV) as a function of QD-laser detun- 
ing (with Acx — 0.5 meV), for different numbers of truncated 
photon states at (a) T = 4 K and (b) T = 20 K: two-photon 
(or more, since this value has converged): black dashed line, 
four-photon: grey solid, lighter line, six-photon (or more): 
red solid, darker line, (c-d) Mean cavity photon number for 
a dot driven system as a function of QD-laser detuning (with 
Acx = 3 meV), for different number of truncated photon 
states with (c) T = 4 K and (d) T = 20 K: one-photon: 
black dashed line, two-photon: grey, lighter solid line, and 
three-photon: red solid, darker line. 



approximation overestimates the power broadening in- 
tensity PL for a negatively detuned cavity. Similar con- 
clusions about the need for multiphonon effects were also 
found in the context of the resonance fluorescence spectra 
of a QD driven cavity- QED system^^ and time-dependent 
excitonic Rabi rotation dynamics However, for a cav- 
ity driven system [Fig.[8fc) and (d)], at low temperatures 
(e.g., 4 K), the weak phonon theory can be accurate. Of 
course, other QD material parameters can yield different 
trends in the role and assessment of multiphonon cou- 
pling. Since the phonon correlation functions that we use 
are well known for an IBM model, the full phonon calcu- 
lation of the intensity PL presents the same level of com- 
putational complexity as the one phonon calculations — 
which is a major advantage of the ME formalism above. 

We also study the influence of quantized multiphoton 
processes in the intensity PL of both an exciton-driven 
and a cavity-driven system. In Figs. [9fa)and[9l^b), using 
cavity excitation {r]c = 30 /ieV), we plot Ix as a func- 
tion of QD-laser detuning for various truncations of the 
photon Hilbert state. We find that it is necessary to in- 
clude up to six-photon processes (e.g., a 13 state model) 
to correctly describe the effect of cavity photons; note 
that including more than six photon processes yields an 
identical result to the six photon calculations, so these 
values have converged on the numerically exact answer. 
While the need for sox photon processes may seem sur- 
prising for the relatively small values of g (i.e., 20 /ieV), 
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the exciton-phonon coupling with phonon scattering can 
be more sensitive to quantum cavity- QED effects. In the 
rest of the paper, we thus compute the intensity PL of 
the cavity driven system in a six-photon truncated basis. 

In Fig. [9fc) and [9l^d) we plot the cavity intensity PL 
(/c) for a QD driven system {r]x = 30 /ieV), as a func- 
tion of exciton-laser detuning for various truncations of 
the photon Hilbert state. Unlike a cavity driven system, 
we now find that a two-photon truncation of the photon 
Hilbert space is enough to obtain a correct (i.e., con- 
verged in the photon basis) description of the cavity PL. 
The fundamental difference between a QD-driven and a 
cavity-driven system is due to the fact that the cavity is 
represented as a quantized harmonic oscillator whereas 
a QD is a two- level system. Hence, for our system pa- 
rameters, the cavity is more easily excited into higher ly- 
ing Jaynes-Cummings ladder states, even though the sys- 
tem exhibits relative strong dissipation (i.e., = 2.5g). 
Neglecting quantum multiphoton processes can therefore 
introduce spurious effects in the PL lineshape. A two- 
photon truncation was also found to be sufficient (and 
necessary) for the study of the fluorescence spectrum of 
a QD driven cavity- QED systemic. 

In summary, the results above highlight the need for 
both multiphonon and multiphoton effects for under- 
standing power-broadened intensity PL, even for low 
temperatures (4 K) and rather small cavity-dot coupling 
rates {g = 20 /ieV, cf. = 50 /ieV). 



C. Power Broadening through Coherent Exciton 
Pumping and Cavity Emission 

Next we study the cavity-emitted PL for different in- 
put powers of an exciton pump. In Figs. [TOT a-b) we plot 
the relative cavity mode intensity (/c) as a function of 
exciton-laser detuning for a QD driven system at two 
temperatures, (a) T = 4 K and (b) T = 20 K, for various 
values of r]x- The cw field drives the QD which is now 
detuned to the right (higher energy) of the cavity mode 
by 3 meV. In Fig. [ToT c) we plot the relative cavity mode 
intensity with only ZPL broadening and set {B) = 1 (i.e., 
no coherent or incoherent coupling effects from phonons). 
This then closely corresponds to an atomiclike power 
broadening model, but with the addition of pure dephas- 
ing of the ZPL — a model that is commonly used to ana- 
lyze semiconductor cavity- QED experiments. The power 
broadened intensity lineshape in the absence of phonon 
coupling can be explained by considering two Lorentzian 
lineshapes centered at the two resonances, the relative 
oscillator strengths of which are qualitatively determined 
by their corresponding broadenings. 

Without phonon coupling, as shown in Fig. [TOT c). an 
increasing cw drive causes power broadening of the QD 
exciton which decreases the oscillator strength relative 
to the cavity mode and also excites the cavity resonance 
more. With phonon scattering processes [see Fig. [TOT a)]. 
the first major difference we notice is the apparent nar- 
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FIG. 10: (Color online) (a-b) Relative cavity mode intensity 
as a function of QD-laser detuning for a dot driven system 
at two different temperatures, (a) T = 4 K and (b) T = 
20 K, and for different values of the cw Rabi frequency {rjx = 
20 /ieV: black dashed line, r]x = 40 /ieV: grey, middle solid 
line, and rix = 60 /xeV: red, upper solid line). The cw field 
drives the QD which is detuned to the right of the cavity 
mode by 3 meV {Acx = -3 meV). Note that (B) = 0.91 
at T = 4 K and (B) = 0.73 at T = 20 K. (c) Normalized 
cavity mode intensity in the presence of only ZPL broadening 
with (B) = 1 at T = 4 K. (d) Plot of the x FWHM of the 
intensity PL at the QD resonance as a function of r/c. The 
orange crosses show the FWHM at T = 20 K, the blue circles 
at T = 4 K, and the red inverted triangles show the intensity 
FWHM with only ZPL broadening and (B) = 1 (no phonons). 
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FIG. 11: (Color online) As in Fig.[TO]but with Acx = 3 meV 
[cf. Acx = -3 meV in Fig.fTO]. 



rowing of the QD resonance compared to Fig. [TOTc) 
which is due to the coherent renormalization of the Rabi 
frequency — this manifests in an effective drive whose 
magnitude decreases with increasing temperature; in ad- 
dition, the phonon interactions reduce the dot-cavity cou- 
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pling through g (B) g. Incoherent phonon couphng 
also introduces significant additional broadening of the 
QD exciton due to the >C(<j+) process (see Table |T]), even- 
tually resulting in a new peak near the phonon cut off 
frequency (i.e., at the peak of the spectral bath function 
for phonons, CJ5 = 1 meV). The mean cavity photon num- 
bers (and thus Ic) increase in the presence of phonons 
due to phonon-assisted processes whose magnitude also 
increases with temperature (see Fig. [ToT c). which also 
has a larger ZPL). Comparing the cavity PL characteris- 
tics in Fig.fToTa-c). we see that electron-phonon coupling 
plays a significant role in determining the intensity PL, 
with features that are not at all explained with simple 
atomiclike MEs. 

In Fig. Mid) we plot the x FWHM of Ic at the QD 
resonance as a function of r]x , which is a typical measure- 
ment in experimental studies^. The orange crosses show 
the intensity at T = 20 K, the blue circles at T = 4 K, 
and the inverted red triangles show the FWHM in the 
absence of any phonon coupling. In spite of significant 
reduction of the effective Rabi frequency due to phonon 
coupling at high temperatures (e.g., (5(4 K)) = 0.91 
{B{20 K)) = 0.77), the x FWHM at T = 20 K is higher 
than the FWHM calculated at T = 4 K, which suggests 
that EID more than compensates for phonon-induced 
renormalization of the Rabi frequency. It is noted from 
Table H] earlier that increasing 7 (and thus also , as 
they have the same Lindblad operators) reduces the mean 
cavity photon number. We also note from Table HVl that 
the primary contribution to the broadening are the two 
scattering terms F^^ and F^^ ; these two processes in- 
crease with temperature and driving strength and in- 
troduce additional broadening which reduces /c. At a 
pump rate of rjx ~ 60 /ieV, we see a more rapid increase 
of the X FWHM with phonon coupling due to stronger 
cw laser-phonon coupling. For pump rates greater than 
rjx ~ 40 /ieV, we also see a more rapid increase of the x 
FWHM with phonon coupling. 

In Figs. fTTT a-b) we plot the relative cavity- mode inten- 
sity PL for a dot driven system, where the exciton reso- 
nance is now detuned to the left of the cavity mode by 
3 meV. We again consider two different temperatures of 
the phonon reservoir, (a) T = 4 K and (b) T = 20 K, for 
various rjx. In Fig. fTTT c) we show the normalized cavity 
mode intensity with only ZPL broadening and {B) = 1. 
The generic features of Fig. [11] can be understood ex- 
actly along the lines of the arguments presented above 
for Fig.[TOl However, as discussed earlier, here we obtain 
a significant reduction in the mean cavity photon number 
(less cavity feeding) since the cavity is now energetically 
higher which requires absorption of phonons — compare 
Fig. [TOT a) and Fig. fTTT a). Moreover, the mean cavity 
photon number is smaller than that in the absence of 
phonons due to the renormalized (reduced) dot-cavity 
coupling which, however, increases with temperature. In 
Fig. fTTTd) we plot the x FWHM of the cavity emission 
as a function of rjx which shows very similar features to 
the X FWHM data presented in Fig. [TOl^d). We remark 



that the approximate x FWHM values are extracted nu- 
merically by fitting with a simple Lorentzian model. As 
it is clear from the plots, the power broadening inten- 
sity PL in the presence of phonons are no longer repre- 
sented by a simple system of coupled Lorentzians, and we 
observe pronounced non-Lorentzian lineshapes and clear 
signatures from the phonon bath spectral function. We 
reiterate the point that there is a substantial discrepancy 
at large pumps between the x FWHM from an effective 
Lindblad solution (EPME) and a full polaron ME solu- 
tion (see Table HV|) . This highlights a breakdown of the 
EPME at large pumps which is not too surprising given 
the rather coarse approximations made in its derivation. 
Finally, we comment that the c FWHM (which can be 
extracted from the same PL data) does not show any 
power broadening, which is possibly also a result of the 
reasonably large QD-cavity detunings. 



D. Power Broadening through Coherent Cavity 
Pumping and Exciton Emission 

In this subsection, we focus on a cavity-excited sys- 
tem, where the emitted excited intensity {Ix) is de- 
tected through the QD emission, e.g., through (non- 
cavity mode) radiation modes. In Figs. [I2fa) and [T2lb) 
we plot the relative exciton intensity, Ix^ as a function 
of QD-laser detuning for a cavity driven system at two 
different temperatures, (a) T = 4 K and (b) T = 20 K, 
for various values of r]c. The cw field drives the cavity 
mode which is now detuned to the right of the QD exci- 
ton by 0.5 meV. In Fig.[T2fc) we show Ix with no phonon 
interactions and (B) = 1 (i.e., no coherent or incoherent 
effects from phonons, apart from pure dephasing). We 
discern that the power broadened intensity PL in the ab- 
sence of phonons is quite distinct from the case with finite 
phonon coupling. On the one hand, we have lost the res- 
onance at the phonon spectral function since there is no 
longer a term that involved incoherent excitation through 
the phonon reservoir (F^^ process). On the other hand, 
we see a clear infiuence from phonon-induced exciton 
cavity feeding. In particular, without phonon coupling 
we observe very little emission at the cavity mode fre- 
quency which further demonstrates that phonons play 
a significant role in the power broadened PL lineshape 
through dot-cavity coupling (via cavity feeding). Also 
note that the broadening of the QD resonance increases 
due to enhanced dot-cavity coupling in the presence of 
phonons — which increases with temperature. However 
at higher temperatures the mean exciton number is de- 
creased which is mainly due to phonon-induced reduction 
in ^ ^ {B)g. 

In Fig. \T2\ d) we plot the c FWHM obtained via exci- 
ton emission (Ix)^ as a function of rjc The c FWHM in 
the presence of phonons at low temperatures (T = 4 K) 
is very similar to the case with only ZPL broadening pri- 
marily because of the absence of EID, as the system is 
now cavity-driven, though g (B) g is still temperature- 
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FIG. 12: (Color online) As in Fig. [TOl but for the exciton 
intensity obtained in the presence of a cavity pump. The cw 
drive strengths are rfc — 10,40,60 /ieV. In (d), the orange 
crosses show the c FWHM at T = 20 K, the blue circles at 
T = 4 K, and the inverted red triangles show the FWHM 
with only ZPL broadening and {B) = 1. The cavity-exciton 
detuning is /\cx = 0.5 meV. 



dependent. However, the c FWHM increases with tem- 
perature even though {B)g reduces. We also note that 
the emission at the QD frequency is suppressed with in- 
creasing temperatures. With further increase in tem- 
peratures (i.e., above 20 K — not shown), the intensity 
PL is dominated by emission at the cavity mode. Sim- 
ilar to the exciton pumped system, for pump rates of 
r]c ~ 40 /ieV, we observe a more rapid increase of the 
c FWHM with phonon coupling. For cavity excitation, 
the effective Lindblad solution shown earlier provides a 
very good match for the c FWHM for all rjc pumps val- 
ues studied (see Table |VI]). This is primarily because the 
approximations made for a cavity driven system are less 
restrictive. For instance, there are only two Lindblad 
terms which is unlike the case of a QD-driven system 
where four Lindblad terms are needed (two of which are 
rjx dependent). 

In Figs. [TSTa) and [TSfb) we again plot the relative 
exciton intensity {Ix) as a function of QD-laser detun- 
ing for a cavity driven system, but now with the cavity 
mode detuned to the left (lower energy) of the QD ex- 
citon by 0.5 meV. In Fig. [TST c) we show the normal- 
ized cavity mode intensity with only ZPL broadening of 
the exciton and {B) = 1. We highlight the significant 
difference in the QD emission at the cavity mode en- 
ergy for the two detunings at low temperature [Fig. [I2fa) 
and Fig. [TST a)] — the cavity mode resonance is more pro- 
nounced with a positive /S^cx- This difference is due to 
the asymmetry in the phonon absorption and emission 
probability; at higher temperatures, however, the line- 
shapes becomes more similar. In Fig. fTST d) we plot the c 
FWHM obtained through exciton emission as a function 
of the cw laser Rabi frequency (r/c), which is very similar 



FIG. 13: (Color online) As in Fig. [HI but for Acx = 
-0.5 meV. 



to that in Fig. HSfd). For these studies, the c FWHM 
is mostly independent of the detuning of the dot-cavity 
system. We also mention that x FWHM (which can be 
extracted from the same PL data) only shows marginal 
power broadening , though it may increase with reduced 
QD-cavity detuning. 



E. Integrated Photoluminescence (IPL) 

Finally, we study the integrated photoluminescence 
(IPL). In Figs.EJa) andHUb) we plot the IPL lineshape 
for a dot-driven system as measured through the cavity, 
for the positively detuned cavity mode (Aca^ = 3 meV). 
The total IPL of the cavity and exciton intensity of a 
dot-driven system plotted in Fig. 14(a) does not sat- 
urate with increasing drives — unlike the IPL over the 
Lorentzian centered at the exciton resonance frequency 
plotted in Fig. 14(b) (i.e., the integrated x FWHM, emit- 
ted via the cavity mode), which shows clear saturation. 
The Lorentzian lineshape at the QD exciton resonance 
is computed by subtracting off the background phonon 
sidebands from the intensity PL. The lack of saturation of 
the IPL of the total cavity and QD intensity is attributed 
to the increasing excitation of the background phonon 
continuum with increasing cw drives which is also en- 
hanced at higher temperatures. These theoretical trends 
are consistent with the experimental results of Ates et 

In Figs.fMfc) andfMfd) we show the IPL of the cavity- 
driven system (with /S^cx = 0.5 meV) as measured via the 
QD exciton emission. Neither the total IPL of the cavity 
and exciton intensity of a cavity-driven system plotted in 
Fig. 14(c), nor the IPL over the Lorentzian centered only 
at the cavity resonance frequency plotted in Fig. 14(d) 
(i.e., the integrated c FWHM, emitted via the QD) show 
any saturation. Furthermore, significant power broad- 
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FIG. 14: (Color online) Integrated PL for a dot-driven system, 
(a) and (b), as measured through cavity emission, and the IPL 
of a cavity-driven system, (c) and (d), as measured via QD 
emission; both as a function of square of the cw drive for two 
different temperatures of the phonon reservoir (T = 4 K: blue 
circles , T = 20 K: orange crosses); also plotted are the IPL 
lineshape in the presence of only ZPL broadening and (B) = 1 
(inverted red triangles); the corresponding ^ = 2 /xeV (at 
4 K). (a) The IPL of the total cavity and QD intensity of a dot- 
driven system as measured through cavity emission; Acx = 
3 meV and the intensity PL are shown in Fig. 1111 (b) IPL 
of the cavity mode intensity over the Lorentzian centered at 
the QD resonance frequency, (c) IPL of the QD exciton and 
cavity intensity of a cavity-driven system as measured through 
the QD; here Acx = 0.5 meV and the intensity lineshapes are 
shown in Fig. 1131 (d) The IPL of the QD mode intensity over 
the Lorentzian centered at the cavity resonance frequency. 



ening is found which increases with temperature of the 
phonon bath because of enhanced phonon-mediated cou- 
pling between the dot and the cavity mode. 



V. CONCLUSIONS 

We have presented a detailed analysis of the intensity 
PL power broadening of a semiconductor cavity- QED 
system under coherent excitation of either the QD exci- 
ton or the cavity mode. In particular, we have included 
the interaction of the QD with an acoustic phonon en- 
vironment at a microscopic level, while also accounting 
for exciton-cavity coupling in the regime of cavity- QED. 
We utilized a time-convolutionless ME approach in the 
polaron frame to study the cavity- QED dynamics. The 
interaction of the phonon reservoir with the QD is non- 
perturbative and is limited mainly by the validity of the 
IBM (and the single exciton picture) at high tempera- 
tures. This approach enables us to treat the coherent 
interaction between the QD, cavity and the cw laser field 
to all orders. Central to this approach is the need to ac- 
count for the internal coupling effects which preserves the 
detail balanced condition on the system density operator. 



Our theory also points out some major flaws and restric- 
tions of atomiclike MEs for modeling coherent excitation 
regimes in semiconductor QD-cavity systems. 

Using material parameters close to those measured in 
related semiconductor experiments, various intensity PL 
lineshapes were studied as a function of the excitation 
pump rate, for different temperatures of the phonon bath. 
We computed the full intensity PL lineshape over a range 
of frequencies and extracted approximate Lorentzianlike 
linewidths of the c/x FWHM of the power broadening 
intensity PL for both the QD-driven and cavity-driven 
system. The interaction of the QD with the phonon reser- 
voir is seen to introduce qualitatively different features in 
the intensity PL, especially at high temperatures, and is 
quite distinct from the case of a typical power broaden- 
ing lineshape of a two-level atom. We found a significant 
narrowing of the QD resonance due to the coherent renor- 
malization of the Rabi frequency resulting in a reduced 
effective drive. The interaction with the phonon reservoir 
also reduced the dot-cavity coupling due to renormaliza- 
tion of ^ ^ (B) g. Phonon coupling further introduces 
additional broadening of the QD exciton due to EID and 
results in highly non-Lorentzian lineshapes for the in- 
tensity PL; via incoherent excitation from the phonon 
bath, clear signatures of the phonon bath spectral func- 
tion appear in the exciton driven PL. For a cavity driven 
system, the mean cavity photon number (which is pro- 
portional to the cavity mode intensity), was also found 
to change in the presence of phonon-assisted processes 
which depended sensitively on dot-cavity detuning and 
the temperature of the phonon bath. 

To help explain the underlying physics of electron- 
phonon scattering in these cavity- QED systems, we de- 
rived an effective phonon ME (EPME) of the Lindblad 
form, which facilitates a very simple numerical solution 
to the full ME and also allows the extraction of various 
phonon-induced scattering rates in a physically meaning- 
fully way. In particular, we identified specific phonon- 
mediated processes which cause EID of the QD exci- 
ton resonance and incoherent exciton pumping. We also 
identified Lindblad superoperators which mediate coher- 
ent interactions between the QD and the cavity mode, 
resulting in exciton-cavity feeding — a phenomenon that 
is becoming more familiar in semiconductor cavity-QED 
(e.g., see Refs. 0118^20-22^24-27^29^40). We also studied 
the regimes of validity of the effective Lindblad solution 
and found that, for relatively large QD-cavity detunings, 
the effective Lindblad solution produces a very good fit 
to the full polaron ME solution. However, discrepancies 
can occur for exciton pumped systems with increasing 
drives, where one requires the full polaron ME to obtain 
a more accurate FWHM of the exciton resonance. Using 
these polaronic ME formalisms, we found that phonon- 
induced EID and incoherent coupling between the QD 
and the cavity is fundamental to obtain a complete pic- 
ture of power broadening in these semiconductor systems. 
In particular, our results demonstrate that the cavity 
emitted QD-driven intensity PL can display quite pro- 
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found signatures of the phonon spectral bath function. 
A cavity driven system also contains clear signatures of 
the phonon bath, containing power broadened features 
that are significantly different to those obtained in atomic 
QED. We have demonstrated that pure dephasing of the 
ZPL, frequently cited as being responsible for phonon- 
scattering effects such as cavity feeding, actually only 
plays a very minor role here; the physics of the phonon 
scattering processes that we identity are manifestly dif- 
ferent to the physics of coupled Lorentzian oscillators, 
either with or without pure dephasing of the ZPL. 

The results presented in this paper should be broad 
interest for interpreting current and future experimental 
data from semiconductor cavity- QED systems. A more 
rigorous analysis would use the full polaron ME and a 
less rigorous approach would use the EPME; the latter 
approach allows for a simple and intuitive picture of the 



underlying physics and can be accurate in certain excita- 
tion regimes. Our power broadening results were chosen 
to highlight new regimes where phonon effects are more 
readily visible and can serve to direct experimental focus 
in those regimes. For this work, we have also focused 
our studies to the domain of cw excitation, but the ex- 
tension to pulsed systems is straightforward and will be 
presented elsewhere. 
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Appendix A: Effective Lindblad master equations with simplified electron-phonon coupling 

In this appendix we provide some technical details into the derivation of our effective Lindblad master equation, 
Eq. (p!5]) . We exemplify for the relatively simple case of a cavity driven system, where Xg = hg{a^a~ + a+a) 
and Xu = ihg{a~^a — d^a~). The case of a dot driven system can also be derived using similar steps but with a 
few additional approximations discussed at the end of this appendix. The integrand inside the phonon integral in 
Eq. (HU), dr^^^g^^{Grn{r)[Xrn.e-'^'^y^^/^Xrne'^^^^^ can be approximated (e.g., for m = ^) as 

follows: 

(Al) 



where we have used 



e 



'^'sys^f^X^e'^'sys^/^ - e-^^o^/^X^e^^o^/^, (A2) 



with Hq = hAxL^^d- -h hAcLCi^o.' This corresponds to approximating H^^^ with Hq in the exponential phase and is 
expected to be only valid when the dot-cavity detuning is large compared to ^. It follows that 

e-^Hor/h^t^-^^H[,r/h ^ e-'^-'d)a-, (A3) 
e-^^o^/^a+ae^^o^/^ = e^^^^^a+a. (A4) 

Consequently, this allows us to write an effective Lindblad ME that can (e.g., see Figs[3E]) reproduce the full polaron 
ME solution over a range of dot-laser and cavity-exciton detunings. The major advantage of this approach (used 
on its own or as a compliment) is that it is simpler than the full polaron ME approach and allows one to extract 
various, physically meaningful, scattering processes associated with electron-phonon interactions. Equation ()Aip can 
subsequently be rewritten as 

Gg{T)[Xg,e-'^'^y^^/^Xge'^'^y^^^^p(t)]^H.c. - h^g^ [G^(r)e-^^^- V+aa^a-p(t) + G^(r)e^^^-'a^(j-a+ap(t)] 

- h^g'^ [Gg{T)e-'^^-'o)cr-p{t)cr"'a + Gg{T)e'^^-^ cr"' ap(t)o) a'] 
+ n'g'' [Gl{T)e-'^-'p(t)d)G-G^a + Gl{T)e'^-' p{t)G^ ao) a'] 

- h^g^ [Gl{T)e-'^^-'a^dp{t)a^a- + Gl{r)e'^^-'d^ p{t)a^d] , (A5) 
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where we have neglected the other terms (such as cr^acr^ap^ cr^apcr^a and so on) as they do not contribute to the 
ME. Equation (|A5p can now be simplified to 



h^g^ Re[Gg{T) 
ih^g^ lm[Gg{r)e 
ih^g'^ lm[Gg{T)e 



a^a cr^ap(t) + p{t)a^a a~^a — 2a~^dp{t)a^ a ) 
]{a~^ aa^ a~ p{t) — p(t)cr^ ad) cf~) 
{d^a~a^dp{t) - p{t)d^ a~ &^ d). (A6) 



Using J^^^o u ^rni^) = {B)'^{e^^^^ — 1), the defined scattering rates become, 



while the frequency-shifts, 



ph 



2(5) VKe 



(5) Vim 



dre 



) _ 



(A7) 



(A8) 



which are used to obtain the contribution of the phonon integral in the Lindblad form, Lph(p), and the effective 

Hamiltonian defined in Eq. ([16]). 

For the case of a dot-driven system, we have an additional complication due to the coherent term driving through 

the exciton-phonon bath. Using similar steps as above, and also neglecting contributions involving cross terms between 

operators (j+ (a") and a+a {d^a~)^ which scale as gr]x^ we again obtain an effective Lindblad form of the ME. We 

evaluate the exponential phase terms involving the full system Hamiltonian by replacing H' with HL 
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